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Abstract 

From consideration of the order-parameter distribution, we propose an observable which makes 
a clear distinction between true and quasi long-range orders in the two-dimensional generalized 
g-state clock model. Measuring this quantity by Monte Carlo simulations for q = 8, we construct a 
phase diagram and identify critical properties across the phase-separation lines among the true long- 
range order, quasi long-range order, and disorder. Our result supports the theoretical prediction 
that there appears a discontinuous order-disorder transition as soon as the two phase-separation 
lines merge. 
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The existence of quasi long-range order (LRO) characterizes the critical behavior of the 
two-dimensional XY model [l[ [2], y| as well as its dual, the solid-on-solid (SOS) model to 
describe the roughning transition on a surface {4], Q, E , 0] . By quasi LRO, we mean that 
the spin-spin correlation function decays algebraically, which implies that the system is not 
magnetically ordered. We will refer to a phase having such characteristics as quasiliquid . 
The lack of true magnetic order for the XY model is attributed to spin-wave excitations, 
which are gapless and thus excited at any finite temperature. On the other hand, the quasi 
LRO is broken by the vortex-pair unbinding at the Kosterlitz-Thouless (KT) transition 
which exhibits an essential singularity. Even though the XY model assumes the continuous 
U(l) symmetry in the spin angle 9, essentially the same nature is observed when the angle 
is discretized into q possible values over = 0,^,..., 2&=!2, as long as q is high enough. 
Such a discrete-spin system is called the g-state clock model if two neighboring spins, which 
have 0i = 2^n,i/q and 9j = 2imj/q with integers n« and rij, respectively, interact via cosine 
potential V(9i — 9j) = —J cos(9i — 9j) with a ferromagnetic coupling constant J > 0. One 
can generalize this interaction with preserved symmetry, V(9) = V(—9) = V{9 + 2n), into 
the form given by the Hamiltonian 
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(1) 



where 9^ is the ith spin angle, and the sum runs over nearest neighbors {J. It recovers the 
g-state clock model at p — 1 and approaches the g-state Potts model in the limit of large 
p jlO]. We denote the system defined by Eq. (pQ) as the generalized g-state clock model. 
Since it has been claimed that this model with p = 1 and q > 8 precisely reproduces the 
KT transition we set q = 8 throughout this work. At the same time, the discreteness 
introduces a finite ga p in the spin-wave excitation, making the true LRO realizable at low 



temperatures 
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14|. These two phase transitions are connected by the duality 



15] . While the 



relation, which is exactly established within the Villain approximation 
appearance of the^ quasi LRO is readily detected by observables such as Binder's fourth- 
order cumulant 



tie quas 



or helicity modulus 
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that of the true LRO has been 



20|. It is, however, rather 



observed by changes in specific heat or magnetization 
hard to locate the transition temperature using these quantities, especially for high p values 
where the quasi LRO exists in a very narrow temperature range. Thus alternative quantities 



are required, for example, like a direct observation of the formation of giant clusters 
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FIG. 1: (Color online) Distributions of the order parameter m on the complex plane, obtained for 
p = 1, q = 8, and L = 8. We start from the high-temperature regime and then cool down the 
system slowly. Here are shown three characteristic distributions, where one finds (a) a disordered 
phase at temperature T = 1.50, (b) a quasiliquid phase at T = 0.70, and (c) an ordered phase at 
T = 0.36, where the temperatures are given in units of J/Jcb- 

In this Rapid Communication, we show that the transition can be well localized by a non- 
local order parameter which is obtained from the average spin direction and which makes a 
clear distinction between the true and quasi LROs. Using this quantity it is shown that the 
quasiliquid phase disappears beyond p ~ 2.8, where the transition becomes discontinuous 
just as for the eight-state Potts model. 

Let us consider the generalized eight-state clock model given above on the L x L square 
lattice with the system size N = L 2 . The complex order parameter of this system is defined 
as 

m = iV- 1 ^^ = \m\e i4> . (2) 



As in Ref. 



221 ]. it is instructive to visualize the distribution of m on the complex plane. The 



distributions in Fig. [T] are obtained by running Monte Carlo simulations with the single- 
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cluster update algorithm 



251 ]. and each panel represents a different phase of the 
eight-state clock model at a different temperature. In the leftmost panel [Fig. [11(a)], we see 
the disordered phase in the high-temperature regime. The order parameter m exhibits a 
two-dimensional Gaussian peak around the origin, which may be regarded as a delta peak 
at | m | = in the thermodynamic limit. Figure d^b) illustrates the quasiliquid phase, where 
the order parameter rotates in the <fi direction with nonzero magnitude. Note that both of 
the distributions in Figs. Ufa) and[T](b) manifest a continuous rotational symmetry, which is 
spontaneously broken at a lower temperature as shown in Fig. [0(c). One finds a true LRO 




FIG. 2: (Color online) Double phase transitions for p = 1 and q = 8. (a) The transition between 
quasiliquid and disordered phases is detected by merging of U m curves, (b) The other transition 
occurs between the ordered and quasiliquid phases, which is detected by m^. (c) Extrapolating 
positions of inflection points according to the KT picture, we get T c \ = 0.417(3) and T C 2 = 0.894(1). 
(d) Checking Eq. (jTj) with rj/2 = 0.031, which best describes the data with a single curve. 

being established so that m indicates well-defined directions selected from the eight-fold 
symmetry. 

A major difference between Figs. QJa) and [0(b) lies in the distributions of \m\. The tran- 
sition between the quasiliquid and disordered phases can be detected by means of Binder's 
fourth-order cumulant, 

U m = 1 - - 1 % , 3 
2(|m| 2 ) 2 1 ; 

where (• • •) represents the thermal average [Fig. 12(a)] . The factor of two in the denominator 
of Eq. ([3]) is based on the fact that (|m,| 4 ) = 2 (|m| 2 ) 2 for such a two-dimensional Gaussian 
distribution as in Fig. [1(a). We should note that U m does not detect the transition between 
the ordered and quasiliquid phases since they differ only in the angular direction on the 
complex plane. Henceforth, we need a quantity capturing the change along 0. In the same 
spirit as U m , one may define a cumulant as 



5(0 4 

U * = 1 - 77T2 ' ( 4 ) 
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where <fi = (2tt) 1 {q<p mod 27r) so that goes to zero when the distribution is uniform with 
respect to <fi. Or we may alternatively have 

m = (cos(g0)) , (5) 

which yields a finite value when <p is frozen but again vanishes when <p is isotropically 
distributed [Fig. [2(b)] . Provided that the system is nearly ordered with large enough q, 
we approximately have <ft ~ from Eq. ([2]) so that « (cos(27rn)) with n = 

N^ 1 J2j n j- By duality, the integer field rij can be mapped to a charge distribution in the 



lattice Coulomb gas [15j and the approximate expression for has been introduced in 
Ref. {j} to monitor the fugacity of charged particles under numerical renormalizat ion-group 
calculations. Since the quasiliquid phase exists between the ordered and disordered phases 
for q = 8, we have two separate transitions at T = T c \ and T C 2, which are clearly detected by 
the above quantities. Note the movements of data points in Figs.[2^a) and [2(b) with different 
system sizes. Since the position of an inflection point, T*, would correspond to where the 
transition occurs in the thermodynamic limit, we may extrapolate them according to the 
KT scenario, 

lnL~ |T,-T c p 1/2 , (6) 
to estimate the critical temperatures both for the upper and lower transitions [Fig. [2(c)]. 

I L 

In addition, regarding (\m\) = L^ v ' 2 m(L } T) around T cl [26], we replace the dependency on 
both of L and T by that on a single variable = m<f,(L, T) so that 

(H) = l-'/ 2 /K). (7) 

In other words, plotting LPl 2 against m^, data from different sizes are expected to fall 
on a single curve if one correctly selectsj?. This provides a way to determine rj even without 



precise knowledge of T cl (see also Ref. [27|). The best fit is found at rj/2 = 0.031(2) as shown 



in Fig. [2(d), while the theoretical value is given as r//2 = 1/32 = 0.03125 at T = T cl 141 ] . 

When altering the potential shape by increasing p in Eq. (CQ), one may well expect that the 
two transitions will eventually transform into a single, discontinuous transition at a certain p 
value as the Potts-model limit is approached. How this happens can be found by numerical 
simulations, and a phase diagram thereby obtained is shown in Fig. EJ^a). It seems that the 
two phase-separation lines merge as p approaches 3.0. A better estimate is obtained by look- 
ing at magnetic susceptibility. Recalling that susceptibility x — N(ksT)^ 1 ((\m\ 2 ) — (\m\) 2 ) 
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FIG. 3: (Color online) (a) Phase diagram of the generalized eight-state clock model, (b) Sus- 
ceptibility as a function of T at p = 2.8, where the dotted lines describe \ ~ ~~ Tc| -1 ' 2 > with 
T c = 0.17184. (c) Ve as a function of T. For each p, system sizes are given by L = 16,32,64, 
and 128 from right to left, (d) Estimation of 77 for p = 2.6 by Eq. ([9]). Seventh-order polynomials 
are used to find the best fit, from which 77 = 0.41(4) is estimated. Inset: T C 2 at p = 2.6 against 
(InL) -2 up to L = 512. It moves slower than predicted by Eq. ([6]) as shown in comparison with 
the straight dotted line. 

corresponds to the sum of correlations, we may argue that its divergence implies long-ranged 
correlations over the system, a key feature of the quasiliquid phase. If p is small enough to 
exhibit the quasiliquid phase, susceptibility indeed diverges over a finite temperature range. 
For p = 2.8, however, we find that data points fall on x ~ |T — T C |~ L2 which has only one 
singular point at T = T c [Fig. EJ^b)]. This T c is also consistent with the results from U m and 
777,0. We therefore conclude that the quasiliquid phase shrinks to a single point at p ~ 2.8. 
Furthermore, the distribution of energy per spin, E, exhibits double peaks for p > 2.8. By 
analogy with U m , we introduce the following quantity: 

1 ((E-(E)) 4 ) 



3((E-(E)) 2 ) 2 ' 



(8) 



Recall that if a scalar variable x has a one-dimensional Gaussian distribution with zero 
mean, one readily finds (x 4 ) = 3 (x 2 ) 2 . Consequently, Ve will vanish when there exists 
a single peak positioned at (E). It will approach a nontrivial value, however, when the 

6 



energy distribution has double peaks on opposite sides of the average value (E). A similar 



attempt to define such a quantity has already been made in Ref. 



26] for characterizing a 



discontinuous transition. Figure [3](c) shows that Ve remains finite at p > 2.8, which signals 



a change to a discontinuous transition 
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The concept of universality suggests that the critical properties will be kept the same in 
the vicinity of p = 1. However, one may ask if the natures of the transitions between the 
ordered and quasiliquid phases and between the quasiliquid and disordered phases depend 
on the value of p. Applying Eq. ([7j) to higher p values, we find that rj = 1/16 cannot be ruled 
out even when p approaches 2. However, the quality of fit severely deteriorates at p higher 
than 2, possibly due to that our magnetization data are easily influenced by the proximity 
of the upper transition. On the other hand, with the same motivation as in Eq. ([7]), we may 



characterize the upper transition by means of the following scaling relation 



27|: 



X = L 2 ^g(U m ). (9) 

This method yields rj = 0.24(1) at p = 1.0 in agreement with the prediction of 1/4 = 0.25 
for the KT transition [l|. It is not very surprising that r] tends to be underestimated here if 



taking into account the logarithmic correction involved in susceptibility [28|, |29j. We observe 
from our numerical data that the criticality deviates from the standard KT type below the 
merging point. If we take p = 2.6, for instance, the best fit is found at 77 = 0.41 and the size 
dependence of the transition temperatures deviates from Eq. [Fig. EJ^d)]. Still, it remains 
to be investigated in detail how the critical behavior begins to change or if the standard 
KT behavior is recovered for even larger lattice sizes (L > 512) in spite of the data collapse 
shown in Fig. [3](d) for lattice sizes up to L = 512. 

In summary, we have proposed a practical quantity to distinguish the true and quasi 
LROs, based on the order-parameter distribution. Using this quantity, we have provided a 
phase diagram on the p — T plane for the generalized eight-state clock model. It has been 
shown that a discontinuous transition appears when the phase-separation lines merge into 
one at p ~ 2.8. We have also checked critical properties along the lines and found changes 
in scaling behaviors before reaching the merging point from numerical calculations up to 
L = 512. 
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